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ABSTRACT 

This  thesis  models  the  motion  of  the  NPS  satellite  ORION 
using  the  three  Euler  rotation  angles.  The  simulation 
program  also  models  aerodynamic  drag  and  gravity  gradient 
torques.  Simulations  are  performed  to  analyze  the  effect  of 
changing  the  inertia  ratio  on  satellite  stability  and 
performance.  The  active  nutation  control  method  was  also 
simulated  and  an  example  of  its  operation  given.  The  amount 
of  time  required  to  reduce  nutation  to  an  acceptable  level 
was  found  to  be  dependent  on  the  initial  nutation  angle, 
spin  rate,  and  thruster  size.  The  slower  the  spin  rate,  the 
greater  the  torque  that  could  be  used  during  a  given  firing 
arc  resulting  in  faster  nutation  control. 
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I.   CONTROL  SYSTEM  OVERVIEW 

This  thesis  is  intended  to  act  as  the  preliminary  ground 
work  for  the  attitude  control  system  of  the  Naval 
Postgraduate  School  satellite  ORION.  The  background  for 
ORION  is  detailed  in  Reference  1.  The  orbital  control 
system  deals  with  the  motion  of  the  center  of  mass  in  an 
orbit  around  the  earth,  while  the  attitude  control  system 
deals  with  the  angular  motion  of  a  satellite  around  its 
center  of  mass.  The  attitude  control  system  must 
incorporate  attitude  sensors,  control  laws,  actuators  or 
thrusters,  the  effect  of  disturbing  forces  and  the  dynamics 
of  the  satellite. 

ORION  is  being  designed  to  be  launched  via  the  extended 
Get-Away-Special  (GAS)  cannister  aboard  the  Space 
Transportation  System  (Space  Shuttle) .  ORION  is  currently 
under  development  and  many  of  the  systems  which  will 
determine  its  final  characteristics  have  not  been  fixed. 
For  the  purpose  of  this  thesis,  ORION  is  considered  to  be  a 
cylindrical  satellite  of  uniform  mass  with  dimensions  of 
0.889  meters  height  and  0.4826  meters  diameter.  It  has  a 
mass  of  113  kilograms.   Figure  1.1  shows  a  diagram  of  ORION. 

The  design  of  a  control  system  is  highly  dependent  on 
the  mission  of  the  satellite.  Reference  2  provides  a 
section  on  the  analysis   of   mission   requirements   and   the 
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possible  control  systems  which  will  fulfill  mission 
requirements.  Tables  1  and  2  show  the  trade-offs  that  need 
to  be  made  when  choosing  a  type  of  control  system  [Ref. 
2:pp.  IV-12,IV-12] . 

ORION  is  being  designed  to  act  as  a  payload  bus  for  a 
variety  of  missions.  General  requirements  are  that  the 
satellite  fly  in  orbits  from  200  kilometers  to  2000 
kilometers.  Attitude  sensor  accuracies  will  not  be  less 
than  2.0  degrees.  The  life  of  the  satellite  will  depend  on 
a  mission  payload  versus  fuel  trade-off.  For  missions  that 
do  not  require  a  high  degree  of  pointing  accuracy,  spin 
stabilization  provides  the  simplest  means  of  attitude 
control.  A  spinning  satellite  will  be  less  susceptible  to 
nutation  if  spun  about  a  maximum  moment  of  inertia  and  will 
be  stable  if  energy  can  be  dissipated.  The  equations 
governing  satellite  motion  are  developed  in  Chapter  III. 

The  purpose  of  the  control  system  is  to  maintain  the 
satellite  at  a  specified  orientation  in  space.  If  the 
satellite  is  perturbed  from  this  orientation  by  external 
factors  such  as  gravity  gradient  torques,  aerodynamic  drag 
or  nutation  then  the  control  system  must  sense  the  change 
and  apply  the  control  torques  necessary  to  reorient  the 
satellite.  It  must  also  be  able  to  reorient  the  satellite 
to  accommodate  any  change  in  the  mission  requirements. 
Figure  1.2  shows  a  diagram  of  the  basic  control  system.  The 
feedback  loop  and   error   determination   process   may   occur 
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TABLE  1.2 
ATTITUDE  CONTROL  METHOD  APPLICATIONS 
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either  on  board  the  satellite  or  via  a  satellite  to  earth 
station  telemetry  link.  On-board  processing  requires 
preprogramming  the  reference  angles  into  the  processor  and 
is  limited  to  those  periods  when  the  control  system  needs  to 
act  only  as  a  regulator  and  any  disturbances  are  small.  A 
telemetry  link  is  required  to  control  the  satellite  during 
major  orientation  periods  and  flight  maneuvers.  It  also 
acts  as  a  back-up  for  the  on  board  processor. 

There  are  three  distinct  phases  of  the  attitude  control 
problem.  They  are:  1)  launch,  consisting  of  activities 
required  to  take  the  satellite  into  a  preliminary  orbit;  2) 
acquisition,  in  which  the  satellite's  preliminary 
orientation  and  maneuvers  are  accomplished;  and  3)  mission 
operations,  in  which  the  mission  requirements  are  carried 
out  [Ref .  3 :p.  12] . 

A .   LAUNCH 

Control  during  launch  is  limited  to  the  launch  vehicle 
and  has  little  relationship  to  the  attitude  control  of  the 
satellite.  Of  primary  concern  is  the  detachment  of  the 
satellite  from  the  launch  vehicle  in  some  preliminary  orbit. 
The  satellite  will  have  the  angular  momentum  of  the  launch 
vehicle  plus  some  velocity  imparted  from  the  separation 
process.  ORION  will  be  carried  into  orbit  by  the  Space 
Shuttle  and  launched  from  a  GAS  cannister.  Reference  1 
details  the  cannister  and  the  launch  process.   The  launch 
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vector  is  along  an  outward  radial  from  Earth's  center  at  4 
feet  per  second. 

B.   ATTITUDE  ACQUISITION 

The  acquisition  phase  commences  after  ORION  is  detached 
from  the  Shuttle.  During  this  phase  the  control  problem 
consists  of  determining  the  attitude  of  the  satellite, 
spinning  the  satellite  up  to  the  desired  spin  rate, 
reorienting  the  satellite  to  align  the  maneuvering  thruster 
with  the  velocity  vector  required  to  place  the  satellite  in 
its  mission  orbit  and  maintaining  the  correct  attitude 
during  the  orbit  transfer.  Once  in  its  final  orbit,  the 
control  system  must  reorient  the  satellite  to  the  desired 
mission  attitude. 

NASA  prohibits  the  firing  of  thrusters  over  5  pounds 
sooner  than  45  minutes  after  separation.  Thrusters  under  5 
pounds  cannot  be  fired  closer  than  200  feet  to  the  shuttle. 
At  a  launch  velocity  of  4  feet  per  second,  this  means  the 
satellite  will  drift  for  50  seconds  before  beginning  spin  up 
for  attitude  orientation.  During  this  drift  time  the 
primary  forces  acting  on  the  satellite  are  aerodynamic  drag 
and  gravity  gradient  torque.  Additional  torques  will  be 
present  if  the  satellite  contacts  the  cannister  side  during 
the  launch  process.  The  environmental  torques  are  developed 
in  Chapter  III. 

The  attitude  of  the  satellite  can  be  predicted  during 
the  drift  period  by  applying  these  torques  to  the  basic 
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satellite  equations  of  motion.  Some  assumptions  must  be 
made  concerning  the  initial  conditions  at  the  time  of 
launch.   These  are: 

(1)  The  satellite  does  not  impact  the  cannister  during 
the  launch  process.  If  this  happens  the  satellite 
will  tumble  and  the  attitude  cannot  be  predicted. 
The  motion  would  have  to  be  observed  by  the  Shuttle 
crew  and  reported. 

(2)  The  satellite  is  aligned  as  shown  in  Figure  1.3  with 
the  z  axis  pointing  towards  the  earth  and  the  x  axis 
in  the  direction  of  flight.  The  initial  values  of  $, 
Y  and  0  are  zero. 

(3)  The  aerodynamic  drag  force  acts  only  against  the 
satellite's  center  of  pressure  and  opposite  to  the 
direction  of  flight. 

After  the  satellite  is  200  feet  from  the  shuttle,  it  is 
spun  up  to  60  rpm  and  the  actual  attitude  determined  by  the 
horizon  and  sun  sensors.  A  description  of  the  horizon  and 
sun  sensors  is  given  in  Chapter  II.  Once  the  orientation  of 
ORION  is  determined,  control  torques  may  be  applied  to 
orient  it  to  any  desired  attitude.  During  this  time  system 
tests  may  be  performed  to  verify  the  correct  operation  of 
the  satellite.  Details  concerning  control  torques  are  found 
in  Chapter  III. 

The  next  portion  of  the  acquisition  phase  is  to  orient 
ORION  for  orbital  maneuvering.  It  is  unlikely  that  the 
Shuttle  will  be  able  to  launch  ORION  into  its  final  orbit. 
For  its  first  mission  ORION  will  be  required  to  extend  four 
magnetometers  from  the  body.  This  deployment  will  have  two 
major  effects.  It  will  despin  the  satellite  and  it  will 
increase   the   moments   of   inertia   about   the   spin   and 
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Figure  1.3   Initial  Axes  Orientation  and  Attitude  Angles 
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transverse  axes.  Once  the  magnetometers  are  deployed,  the 
satellite  will  need  to  be  oriented  to  align  the  propulsion 
thruster  with  the  velocity  vector  which  will  maneuver  ORION 
into  its  final  orbit.  Details  on  the  final  orbit  are  not 
yet  available. 

C.  MISSION 

The  third  phase  of  the  attitude  problem  is  the  control 
of  the  satellite  while  it  is  performing  its  mission. 
Generally  this  reguires  the  control  system  to  act  as  a 
regulator  to  maintain  the  correct  orientation.  If  the 
reguired  attitude  is  fixed,  the  control  loop  may  be 
programmed  into  a  micro-processor  and  function  autonomously. 
The  telemetry  link  would  be  used  to  monitor  the  processor 
and  act  as  a  back-up.  Control  by  the  earth  station  may  be 
necessary  to  orient  the  satellite  for  maneuvers  or  to 
correct  for  a  large  disturbance.  Figure  1.4  summarizes  the 
control  problem  in  algorithmic  form. 

D.  ASSUMPTIONS 

The  purpose  of  this  thesis  is  to  form  a  model  of  ORION 
and  simulate  its  motion  during  the  three  phases  of  the 
control  problem.  It  will  model  the  satellite's  eguations  of 
motion,  the  environmental  disturbance  torgues,  and  the 
control  torgues  and  study  the  interactions  between  these 
factors.  As  stated  earlier,  ORION  is  still  in  the  early 
stages  of  design.   Many  of  the  systems  which  determine   the 
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flight  characteristics  have  not  been  developed.  It  is 
necessary  therefore  to  make  a  number  of  assumptions  which 
simplify  the  simulation  problem.  As  the  design  of  ORION 
matures  assumptions  may  be  removed  and  the  model  updated. 
This  will  help  in  determining  design  trade-offs  and  increase 
the  accuracy  of  the  model. 

The  first  major  assumption  is  that  the  satellite  has  a 
uniformly  distributed  mass.  This  locates  the  center  of  mass 
at  the  origin  of  the  body  aligned  axis  as  shown  in  Figure 
1.1.  The  moments  of  inertia  are  dependent  on  the  deployment 
of  the  magnetometers,  but  in  the  launch  configuration  Ia  is 
3.239  kg-m2  and  I-j-  is  8.916  kg-m2 .  As  the  mass  distribution 
changes  because  of  component  location  or  changes  in  fuel 
weight,  payload,  etc.,  these  moments  of  inertia  will  change. 
If  the  moments  of  inertia  are  located  along  the  body  aligned 
axes,  the  cross  products  of  inertia  are  zero.  The  second 
major  assumption  is  that  the  orbit  will  be  circular.  This 
results  in  a  constant  angular  velocity  w0  around  the  earth. 
The  radius  of  orbit  also  remains  constant  throughout  the 
orbit. 

The  model  of  the  satellite  will  first  be  constructed 
assuming  the  satellite  is  a  rigid  body  with  no  internal 
energy  effects  such  as  sloshing  fuel  tanks.  It  will  also  be 
assumed  that  there  will  be  no  significant  changes  in  the 
moments  of  inertia  because  of  the  control  thruster  burns. 
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Energy  effects  can  be  taken  into  account  later  in  the 
simulation  process. 

The  model  will  also  assume  that  the  spin  rate  of  the 
satellite  will  remain  constant.  In  actuality  the  drag  force 
on  the  satellite  will  cause  the  spin  rate  to  decrease 
significantly  in  approximately  50  days  [Ref.  4:p.  34].  This 
will  be  ignored  in  the  model  since  if  a  significant 
reduction  is  noticed,  then  the  spin  thrusters  can  easily  be 
fired  to  spin  the  satellite  without  affecting  the 
orientation. 

When  a  satellite  is  in  orbit,  it  is  subjected  to  a 
variety  of  environmental  forces  and  moments.  These 
environmental  disturbances  include  gravity  gradient  torque, 
aerodynamic  drag,  solar  radiation  pressure,  magnetic  field 
torques  and  micrometeor  and  radiation  impacts.  Gravity 
gradient  and  aerodynamic  drag  effects  are  dependent  on  the 
radius  of  the  orbit.  As  the  radius  increases  these  effects 
are  decreased.  Solar  pressure  and  magnetic  torques  are 
important  over  a  long  period  of  time,  but  can  be  ignored 
when  looking  at  the  short  term  motions.  The  primary  effect 
of  these  environmental  disturbances  is  noticed  in  the 
nutation  and  precession  of  the  spinning  satellite  over  time. 
This  model  will  assume  that  the  total  torque  on  the 
satellite  will  be  the  sum  of  gravity  gradient  torque, 
aerodynamic  drag  torque  and  applied  control  torques.   The 
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model  will  allow  any  additional  torques  to  be  applied  at  the 
simulator's  discretion. 
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II.   ATTITUDE  SENSORS 

Before  control  can  be  applied  to  a  satellite,  the  actual 
orientation  in  space  must  be  determined.  This  is  done  by 
the  attitude  sensor  package.  Many  devices  and  methods  can 
be  used  as  attitude  sensors.  These  include  rate  gyroscopes, 
star  sensors,  sun  sensors,  earth  limb  or  horizon  sensors, 
and  magnetometers.  This  chapter  will  describe  how 
magnetometers,  sun  sensors  and  horizon  sensors  operate. 
These  three  types  of  sensors  are  currently  planned  for  use 
in  ORION. 

A.   SUN  SENSORS 

Sun  sensors  are  used  to  measure  the  angle  of  the  sun 
relative  to  the  spin  axis  of  the  satellite.  The  sensor 
consists  of  a  pattern  of  photoelectric  cells  and  a  housing 
that  limits  the  angle  at  which  the  sun's  rays  contact  the 
cells.  Figure  2.1  shows  a  schematic  of  the  sun  sensor  [Ref. 
3:p.  163]. 

The  field  of  view  slit  is  oriented  parallel  to  the  x-z 
plane  and  perpendicular  to  the  spin  axis.  The  angle  3  de- 
termines the  amount  of  sunlight  on  the  photo  cells  which  in 
turn  produces  a  given  amount  of  voltage.  As  the  spin  axis 
moves  in  inertial  space,  usually  due  to  a  torque,  3  changes 
and  the  voltage  from  the  photo  cells  changes.  The  accuracy 
of  the  sun  sensor  is  affected  by  the  quality  of  solar  cells 
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used,  the  presence  of  any  highly  reflective  components  of 
the  satellite  and  the  location  of  the  sensor  on  the  body.  A 
typical  accuracy  is  0.01  degrees  [Ref.  3: p.  17].  The  rate 
of  change  of  3  is  determined  by  calculating  the  time  rate  of 
change  of  the  output  voltage. 

Sun  sensor  data  can  be  used  to  calculate  nutation  and  to 
calculate  the  yaw  angle  H1 .  The  method  for  determining 
nutation  is  described  later  in  this  chapter,  4*  is  deter- 
mined by  using  information  from  the  sun  sensor  and  the 
horizon  sensors.  The  motion  of  the  spin  axis  is  a  combina- 
tion of  both  $  and  ¥.  If  $  is  known  from  the  horizon 
sensors  then  for  small  angles, 

^2   =   62  _  $2  (2.1) 

and 

•      •    • 

V   =   6  -  $  (2.2) 

B.   HORIZON  SENSORS 

Earth  limb  or  horizon  sensors  are  IR  sensors  with  a  set 
field  of  view  that  measures  the  difference  in  temperature 
between  space  and  the  earth's  limb.  Space  is  considered  to 
have  a  temperature  of  4  degrees  Kelvin  and  the  Earth  has  a 
temperature  of  300  degrees  Kelvin.  As  the  sensor's  field  of 
view  sweeps  from  space  onto  the  earth  a  pulse  is  started. 
The  pulse  ends  as  the  FOV  goes  from  earth  back  onto  space. 
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Two  sensors  are  located  equidistant  from  the  center  line 
of  the  satellite  as  shown  in  Figure  2.2.  As  the  satellite 
spins,  the  two  sensors  create  a  pair  of  pulses  as  they  cross 
the  earth.   The  pulse  width  is  determined  by 

PW   =   R/2h6  (2.3) 

where  R  is  the  earth  chord  length  (Figure  2.3),  h  is  the 
altitude  of  the  satellite  and  0  is  the  spin  rate  of  the 
satellite. 

The  sensors  can  determine  $,  <3>  ,  and  9.  The  spin  rate  0 
can  be  computed  by  measuring  the  time  between  successive 
pulse  starts.  $  is  determined  by  comparing  the  pulse  widths 
of  sensor  1  and  sensor  2.  When  the  pulse  width  of  sensor  1 
is  less  than  the  pulse  width  of  sensor  2,  $  is  positive. 
When  the  pulse  width  of  sensor  1  is  greater  than  that  of 
sensor  2,  $  is  negative.  This  relationship  is  shown  in 
Figure  2.3.  The  value  of  $  is  computed  by  comparing  the 
magnitudes  of  the  pulse  widths  with  those  of  $  =  0  and  using 
the  methods  in  Chapter  11  of  Reference  5.  $  is  determined 
by  comparing  the  pulse  widths  of  successive  pulses  over 
time. 

The  accuracy  of  the  sensor  is  dependent  on  the  quality 
of  components,  and  the  presence  of  infrared  sources  other 
than  the  earth.  Masks  are  used  to  block  the  effects  of  the 
sun  and  avoid  the  saturation  of  the  sensors.  The  satellite 
itself  may  act  as   a  noise   source,  however  shielding  and 
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Figure  2.2   Horizon  Sensor  Location 
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Figure    2.3      Relationship   Between   Earth   Chord   and  Roll   Angle 
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knowing  the  infrared  signature  can  reduce  the  effects. 
Pages  365-409  of  Reference  5  give  computational  methods  used 
to  eliminate  errors  and  determine  the  accuracy  of  the 
sensors.  In  general  the  accuracy  of  a  horizon  sensor  can  be 
assumed  to  be  0.1  degrees  [Ref.  3:p.  17]. 

C.   MAGNETOMETERS 

Considering  the  accuracy  in  attitude  measurements 
reguired  for  ORION,  the  sun  and  horizon  sensors  will  suffice 
for  the  sensor  package.  Since  magnetometers  are  going  to  be 
used  in  the  mission,  they  can  be  used  to  provide  additional 
attitude  inputs.  These  inputs  can  be  compared  with  those 
obtained  from  the  sun  and  horizon  sensors  to  help  eliminate 
noise  and  errors  in  measurements. 

For  attitude  sensing  purposes,  three  mutually  orthogonal 
magnetometers  are  reguired.  These  may  either  be  internally 
mounted  or  deployed  on  booms  [Ref.  3:p.  250].  The 
magnetometers  consist  of  a  coil  around  a  ferromagnetic  core 
of  known  permeability  and  area.  When  in  the  presence  of  a 
magnetic  field,  a  voltage  is  produced  that  is  proportional 
to  the  field  component  along  the  solenoid  axis  and  is  given 
by 

V   =   -ANy(dBj/dt)  (2.4) 

where  B^  is  the  field  component  [Ref.  3:p.  181].  By 
combining  the  field  components  from  the  three  magnetometers, 
the  solenoid  axis  directions  can  be  computed.    If  the 
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solenoid  axes  are  aligned  with  the  body  axes,  the  attitude 
angles  can  easily  be  found. 

Errors  in  the  magnetometers  may  be  the  result  of  the 
model  used  to  predict  the  earth's  magnetic  field,  the 
influence  of  the  satellite's  electrical  components, 
crosstalk  between  the  magnetometers  or  a  misalignment  of  the 
solenoid  axis.  At  high  orbital  heights,  the  exact  magnitude 
and  direction  of  the  earth's  magnetic  field  is  not  well 
known  and  the  model  used  to  predict  the  magnitude  and 
direction  of  the  flux  lines  may  include  considerable  errors. 
At  altitudes  above  1000  km,  the  field  strength  decreases  as 
the  inverse  cube  of  the  distance  and  the  magnetic  field  of 
the  satellite  may  dominate.  Crosstalk  involves  the 
influence  of  one  magnetometer  on  the  others  as  a  result  of 
eddy  currents  that  are  produced.  Misalignment  may  occur  if 
the  solenoid  axis  is  not  aligned  with  the  body  axis.  This 
may  occur  upon  deployment  or  if  there  is  some  deformation  of 
the  boom  during  a  spinning  maneuver.  [Ref.  3:  pp. 
181,182,250] 

D.   NUTATION 

All  sensor  measurements  described  so  far  have  assumed 
that  the  satellite  was  not  nutating.  If  nutation  occurs, 
the  sensor  data  will  exhibit  an  oscillatory  behavior.  The 
frequency  of  oscillation  is  the  same  as  the  body  nutation 
rate  0Jn.   This  motion  can  be  used  to  calculate  the  nutation 
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angle  directly  from  sensor  data.  The  derivation  and  details 
are  given  on  pages  539-548  of  Reference  5. 

If  the  nutation  is  monitored  by  a  sun  sensor,  then  the 
maximum  nutation  angle  is  given  by 

n   =   63/2R3  (2.6) 

where  <53  is  the  range  of  nominal  sun  vectors  and  R  is  the 
ratio  of  observed  sun  angle  variations  to  the  maximum 
nutation  amplitude.  R  is  dependent  on  the  satellite's 
moments  of  inertia  and  the  location  of  the  sun  sensor  [Ref. 
3:p.  250]. 

Knowing  the  nutation  angle  and  the  attitude  angles,  the 
satellite  can  determine  when  to  fire  control  thrusters  to 
reduce  nutation  and  reorient  itself.  This  can  either  be 
done  autonomously  or  via  the  telemetry  link.  Examples  of 
active  nutation  control  and  orientation  are  given  in  Chapter 
V. 
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III.   ATTITUDE  DYNAMICS  AND  CONTROL 

A.   ATTITUDE  DYNAMICS 

Newton  and  Euler  first  described  the  motion  of  a  system 
of  particles  in  inertial  and  rotational  frames  of  reference. 
These  are  described  and  detailed  in  a  number  of 
publications.  This  chapter  will  present  a  summary  of  these 
derivations  as  they  relate  to  modelling  ORION' s  motion. 

The  satellite  operates  with  respect  to  a  reference  frame 
which  rotates  around  the  earth  in  the  orbital  plane.  As  was 
shown  in  Figure  1.3,  the  x  axis  is  in  the  direction  of 
flight,  the  z  axis  points  towards  the  earth  and  the  y  axis 
is  perpendicular  to  the  orbital  plane. 

The  orbital  axes  are  related  to  the  fixed  reference 
frame  of  the  earth  by  the  direction  cosine  matrix  C.  C  is  a 
linear  combination  of  Euler  rotations  defined  as 

where  Fb  is  the  body  reference  frame  and  F ±    is  the  inertial 
reference  frame  [Ref.  5:p.  23].   The  three  rotations  are 


C(<|>)       = 


10  0 

0  cos  <p        sin  $ 

0        -sin  <J>        cos  <$> 


(3.1) 
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C(0)   = 


C(V)   = 


COS  0 

0   -sin  0 

0 

1       0 

sin  0 

0     cos  0 

cos  ¥ 

sin  ¥    0 

-sin  m 

cos  W    0 

0 

0     1 

(3.2) 


(3.3) 


The  most  commonly  used  sequence  is  the  3-2-1  sequence  which 
results  in 


C   = 


cosfcosg  sin^cos^  sihFsincf) 

-siny  cos0sincj)  +cos^  sin0  cos<|) 

-sinycoso  cosTcosc})  cos^sincf) 

+siny  sinGsincj)  -cosjj  sin¥  sin0 


-sine 


-sind)Cos0 


COS0COS( 


(3.4) 


The  satellite  orientation  is  related  to  the  orbital  axes 
by  the  body  fixed  coordinate  axes  (x,y,z).  The  deviations 
of  the  body  axes  and  the  orbital  axes  are  qiven  by  the 
angles  tj>  (phi)  ,  ¥  (psi)  ,  and  0  (theta)  . 

There  are  six  Euler  equations  that  govern  satellite 
motion  [Ref.  5:p.  32].   They  are: 
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p 

= 

03  *p  + 

f 

p 

= 

mv 

■c 

= 

V 

h 

= 

-oo*hc 

+  T 

h 

= 

1(0 

C 

— 

-w*C 

(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 
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where  co  *    is   defined   as 


co 


*      = 


co 


-co 


'CO, 


CO 


CO, 


-co 


0 


(3.11) 


Equations  (3.5),  (3.6),  and  (3.7)  are  the  translation 
equations  of  motion  and  are  not  discussed  in  this  thesis. 
Equations  (3.8),  (3.9),  and  (3.10)  apply  to  the  attitude 
control  of  ORION. 

The  moment  of  inertia  about  an  axis  is  the  sum  of  the 
masses  times  the  square  of  the  distance  of  the  mass  from  the 
axis.   The  moment  of  inertia  matrix  I  is  given  as 


I   = 


■yx 


zx 


•xy 


-yz 


xz 


■yz 


(3.12) 


where  Ix  =  /  x2  dm,  Iy  =  /y2  dm,  Iz  =  /z2  dm.  Since  ORION 
is  a  symmetrical  body  with  the  body  axes  centered  at  the 
center  of  mass,  the  cross-products  of  inertia,  Ixy,  *yz  and 
Izx  are  zero. 

Expanding  Equation  (3.8)   into  x,  y,  and  z  components 
yields 


hx  =  wzhy  "  wyhz  +  Tx 
hy  =   wxhx  -  wzhx  +  Ty 

hz   =   wyhx  -  wxhy  +  Tz 


(3.13) 
(3.14) 
(3.15) 
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Expanding   Equation   (3.9)   and   substituting   into   (3.13), 
(3.14),  and  (3.15)  and  solving  for  T  yields 


Tx  =  Ixajx  +  (Iz~Iy)a)ycoz  (3.16) 

Ty    =     IyOOy   +   (Ix_Iz)WXWZ  (3.17) 

Tz   =   Izojz  +  (Iy-Ix)wywx  (3.18) 


Tjl  is  the  sum  of  the  external  torques  on  the  satellite  in 
the  ith  axis.  These  torques  include  environmental  torques, 
control  torques  and  any  other  torques  which  may  affect  the 
satellite.   Torques  will  be  discussed  later  in  the  chapter. 

Since  ORION  is  symmetric  about  the  x  and  z  axes,  Ix  =  Iz 
=  1-t,  the  transverse  moment  of  inertia.   The  y  axis  is  also 
the  spin  axis  and  Iy  =  Ia.   The  angular  velocity  vector  co 
decomposes  into  its  Euler  angle  components:   [Ref.  4:p.  30] 

(3.19) 
(3.20) 
(3.21) 

ORION  will  spin  such  that  the  spin  axis  is  in  the  negative  y 
direction.   This  means  that 

ay  =  co0  -  0  (3.22) 

In  the  absence  of  external  torques  the  motion  equations 
reduce  to 

Ita)x  =   (*a~*t)  toy^z  (3.23) 

Itwy   =0  (3.24) 
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wx 

=    4> 

-    ^o 

U)y 

=     0 

"    wo 

WZ 

=    V 

+       (f)(j0o 

ItaJz   =   (It"Ia)a)ya)x  (3.25) 

Equation  (3.24)  implies  that  oo  is  a  constant.  Differen- 
tiating Equations  (3.23)  and  (3.25)  and  substituting  f  or  cj  x 
and  Uy  yields 

(Ia"It) (Ia"It)   2 

oo  H ■ ~ coco=0  (3.26) 

x  2        y  x 

t 

(I  "I.  )  (I  "I.  )    n 

o,  +   a  fc  -  a  fc  (Ao   =  0  (3.27) 

z         2       y  z 

t 
For  stability,  the  coefficients  of  oox  and  ooz  must  be 
positive.  The  y  axis  is  the  maximum  moment  of  inertia  when 
Ia  >  1-t  and  is  the  minimum  moment  of  inertia  when  I-j-  >  Ia. 
If  there  is  no  energy  dissipation,  the  satellite  is  stable 
if  spun  about  either  the  maximum  or  minimum  moment  of 
inertia.  If  there  are  external  torques  or  energy 
dissipation  means,  the  satellite  is  stable  only  if  spun 
about  the  maximum  moment  of  inertia  [Ref.  4:p.  112]. 

Since  there  are  external  torques  present,  these  must  be 
taken  into  account  when  calculating  the  attitude  of  the 
satellite.  Substituting  for  oj  in  Equations  (3.16),  (3.17), 
and  (3.18) , 

Tx  =  It(<t>-0)o¥)  +  (It-Ia)  (wo"0)  (^+^) 

=      It<})-Itco0y  +  (It-Ia)  (    ooo^+ooo^-GY-oOoOt})) 

=      Iti+(It-Ia)  (u3o-u)o0H+((It-Ia)  (w0^)-It^o)^  (3.28) 
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Ty   =   -Iy 0  (3.29) 

Tz   =   It(Y  +  co0J)  +  (Ia-It)  (03O-6)  (i  -coof) 

**  •  •     O      •  •       • 

=    Itvi/+(Ia-It)  (co0e-coo)  necia-^^t)  (ooO-0)  +Itwo)  i 

=   It^+(It-Ia)  (Wo-^oQ),l/-((It-Ia)  (ajO-0)-It^o)i    (3.30) 

Solving  for  the  second  derivatives  of  the  Euler  angles 
yields 

*  =   (Vdt^a)  (wo-wo0H-((It-Ia)  (coo-0)-Itcoo)Y)/It  (3.31) 

Y   =   (^-(Ifla)  (coo-ooo0)^+((lt-la)  (aJo-0)-Ita)o)(}))/It  (3.32) 

0   =   -Ty/Ia  (3.33) 

If  co0  and  0  are  assumed  to  be  constant  the  equations  can  be 
written 

<}>=TX-A<j>-B1'  (3.34) 

4f  =   Tz  -  Ay  +  B  (J)  (3.35) 

0  =   -Ty/Ia  (3.36) 

where: 

A  =   (1  -  Ia/It)  (ooo-0)coo 

B   =   (1  -  Ia/It)  (coo"0)  ~  w0 

The  block  diagram  form  of  the  equations  is  given  in  Figure 
3.1.  As  can  be  seen  from  Figure  3.1,  $  and  4*  are  cross- 
coupled  through  their  rotation  rates. 
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Figure  3.1   System  Block  Diagram 
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The  total  angular  momentum  vector  H  can  be  expressed  as 

~H      =      Ia^y  +  It^t  (3.37) 

2     2      2        . 
where  c^  =  oox  +  u)z  and  lies  in  the  x-z  plane.   oo^  rotates 

relative  to  the  satellite  with  the  angular  velocity  A  where 

*  =   (Ia/It  "  !)wy  (3.38) 

The  body  nutation  frequency  is  defined  as 

w n   =   Ia/It  wy  (3.39) 

The  motion  of  the  satellite  consists  of  the  satellite 
rotating  about  its  spin  axis  and  the  spin  axis  rotating 
about  H  with  the  angular  velocity  ton  [Ref.  4:p.  115].  The 
nutation  angle  n  is  defined  as 

-1  Xtwt 
n  =  sin   (-^p)  (3.40) 

Figure  3.2  shows  the  relationship  between  w^-,  con,  ajy  and  H. 

The  inherent  stability  of  the  satellite  can  be  found  by 
analyzing  the  transfer  functions  of  the  satellite.  Setting 
Tx  and  Tz  equal  to  zero  and  dividing  by  1-j.  yields 

•  •  • 

4>+(l-Ia/It)  (wo-0)GOd})+((l-Ia/It)  (u0-G)-(i)0)V   =   0      (3.41) 
Y+(1-Ia/It)  (u>o-0)a>J-((l-Ia/It)  (Wo-Q)^oH   -   0      (3.42) 

Noting  that  -A  =  (1-Ia/It)  (co0-0)  and  taking  the  Laplace 
transforms  of  the  equations  yields 

c})(s)s2  -  Aco0<i)(s)  "  (A+co0)s^(s)   =   0  (3.43) 
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Figure    3 . 2      Nutation   Cone 


43 


¥(s)s2  -  Xco0^(s)  +  (X+wo)S(f)(s)   =   0         (3.44) 

Solving  for  ^(s)  and  substituting 

,2  2 
„         (A+w  )  s 

(j)(s)  (s  -Aojo) y- <j)(s)   =   0  (3.45) 

(s  -kQ) 

(J>(s)  (s2-  XcuQ)2    +    (  A+coQ)2s2  <|>(s)       =      0  (3.46) 

4>(S)  (s4+(  (X  +  co0)2-Aco0)s2    +    A2co0)       =       0  (3.47) 

Partitioning 


,2  2 

A    d)  ,  _       .  _.  . 

,,    >  -,     , o (3.48) 

<t>(s)      =     1  +  -5 — 9 o 

s    (sZ+((A+ojo)    -Acoq)) 

The  open  loop  transfer  function  is 


^2  2 
o (3.49) 

s2(s2  +((A+o)q)2  -Xa»o) 
The  transfer  function  has  two  poles  at  the  origin  and  two 


poles   on   the   imaginary   axis   at   +j  <J  ( A  +  co0 )  2 -  Au)0   and 
~ J \A A +  wo )   Acoo •   Therefore,  the  stability  of  the  satellite 
depends  on  spin  rate,   inertia  ratio  and  orbit  altitude. 
Table  3.1  shows  the  root  locations  for  different  a  and  coy. 

Because  there  is  no  cubed  term  in  the  transfer  function, 
the  roots  will  always  be  on  the  imaginary  axis.  This 
implies  that  the  satellite  will  tend  to  reach  a  new  steady- 
state  value  for  $  and  ¥  if  perturbed,  rather  than  return  to 
an  initial  state. 
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TABLE  3.1 

ROOTS  OF  THE  TRANSFER  FUNCTION 
(co0  =  -001184) 


Root  Location 

COy 

0 

X 

on  Imaginary  Axis 

6.283 

0.3635 

3.999 

4.00 

1.4882 

-3.067 

3.066 

2.094 

0.3635 

1.333 

1.334 

1.4882 

-1.022 

1.022 

1.047 

0.3635 

.6664 

.667 

1.4882 

-  .5111 

.5106 

B.   ENVIRONMENTAL  TORQUES 

Before  the  motion  of  ORION  can  be  accurately  modelled, 
the  possible  external  disturbances  must  be  studied.  The 
major  disturbances  experienced  by  a  satellite  in  low  earth 
orbit  are  aerodynamic  drag  and  gravity  gradient  torques. 

Aerodynamic  drag  is  a  result  of  the  motion  of  the 
satellite  through  the  upper  atmosphere.  It  can  be  expressed 
as 


1     o 
D   =  y  P  V2  CD  A  (3.50) 


where: 


p  =  density  of  air  at  orbit  altitude 

V  =  velocity  of  satellite  along  orbit 

Cq  =  drag  coefficient 

A  =  surface  area  over  which  the  drag  acts. 
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where: 


The  drag  on  the  body  can  be  represented  by  a  single  force 
acting  on  the  body's  center  of  pressure.  If  this  center  of 
pressure  is  located  away  from  the  center  of  mass,  a  torque 
is  developed. 

The  velocity  of  the  satellite  is 

V   =   wqr  (3.51) 

wo      =      ^e/R3  (3.52) 

V2       =       ye/R3*R2       =        Vfe/R  (3.53) 

ye   =  gravitational  constant 

R   =  distance  of  satellite  from  earth's  center. 

The  density  of  air  decreases  as  the  altitude  increases.  At 
a  200  km  orbit  the  air  density  is  approximately  7.0  x10"10 
kg/m3 .  The  surface  area  is  the  cross-sectional  area  of  the 
satellite;  for  ORION,  A  =  0.429  m2 .  Drag  coefficients  vary 
with  shape,  material  smoothness,  etc.  A  typical  value  for 
this  type  of  satellite  is  Cq  =  2.  The  drag  can  be  expressed 
as 


=  1.71  xlQ14  p  (3>54) 

R 


The  center  of  pressure  for  ORION  is  assumed  to  be 

located  0.076  m  along  the  positive  y  axis.  Since  the  drag 

acts  opposite  to  the  direction  of  flight,  the  total  drag 
torque  on  ORION  can  be  expressed  as 
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TD   =   D  *  L  *  Cos  ^  (3.55) 

where  L  is  0.076  m. 

The  second  major  disturbance  is  due  to  the  gravity 
gradient  that  exists  as  the  distance  from  the  earth 
increases.   The  gravity  gradient  torque  is  defined  as 

G   =  — 5  /(r  xr)  (r.R)  dm  (3.56) 

R 
o 

where: 

R  =   position  vector  of  the  satellite  with  respect 
to  the  earth;  and 

r   =   coordinate  of  an  elemental  mass  from  the 
center  of  mass.   [Ref.  4:p.  131] 

Using  the  3-2-1  Euler  rotation  between  the  body  axes  and  the 
orbit  fixed  axis,  the  gravity  gradient  torques  are 

2 
Gx   =   3co0[  (It-Ia)  sine})   cos  <J>  cos2©]         (3.57) 

Gy   =   0  (3.58) 

2 
Gz   =   3oj0[  (It-Ia)  sinO  cos  6  sin  <J>  ]         (3.59) 

Thus  gravity  gradient  torques  depend  on  orbital  altitude, 
moments  of  inertia  and  satellite  orientation. 

The  equations  describing  aerodynamic  drag  and  gravity 
gradient  torques  can  now  be  used  in  simulating  their  effect 
on  satellite  motion.  The  analysis  of  this  effect  is 
contained  in  Chapter  IV. 
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C.   COMPUTER  MODELS 

The  computer  model  was  written  in  the  Dynamic  Simulation 
Language  (DSL) .  The  program  listing  is  contained  in  the 
Appendix.  The  program  was  intended  to  be  as  general  as 
possible  to  allow  it  to  apply  to  a  variety  of  parameter 
changes.  Of  special  interest  are  I^-,  Ia,  orbital  altitude, 
air  density,  and  spin  rate. 

The  three  Euler  equations  are  integrated  once  to  obtain 
(J)  ,  y  and  9  .  These  are  integrated  again  to  obtain  cj),  ^  ,  and 
0 .  The  initial  values  depend  on  the  situation  being 
analyzed.  Once  the  attitude  angles  have  been  calculated, 
they  can  be  used  to  calculate  wt,  w  n,  \  and  the  nutation 
angle. 

Also  included  is  the  equation  governing  magnetometer 
deployment.  It  is  implemented  in  the  form  of  a  ramp 
function.  The  starting  time  is  the  time  when  the 
magnetometer  is  to  be  deployed.  The  ending  time  is 
calculated  by  dividing  the  desired  boom  length  by  0.0169. 
If  no  boom  length  change  is  desired  the  ramp  function  can  be 
commented  out  of  the  program  and  the  variable  Lj-,  set  to  the 
desired  boom  length.  An  Ly  of  zero  indicates  that  the 
magnetometers  are  in  the  storage  position. 

Changes  in  spin  rates  can  be  accomplished  in  two  ways. 
One  is  to  change  the  initial  condition  of  6  in  a  PARAM 
statement.  This  is  useful  when  evaluating  the  effect  of 
configuration  changes  at  different  spin  rates.   The  second 
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method  simulates  the  use  of  the  spin-up  thrusters.  Ty  can 
be  a  step,  ramp,  or  pulse  input,  or  any  combination  that 
approximates  the  actual  thruster  firing  profile. 

In  addition  to  the  variables  describing  moments  of 
inertia,  spin  rate,  and  thruster  torgue,  many  other 
components  of  the  model  are  left  as  variable  guantities  in 
the  program.  These  can  be  changed  by  PARAM  statements  in 
multiple  runs  to  provide  the  widest  range  of  possible 
parameter  studies. 
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IV.   PARAMETER  STUDIES  AND  ACQUISITION  SIMULATION 

A.  SIMULATION  MODEL 

The  design  of  any  modern  system  involves  the 
investigation  and  analysis  of  how  the  system  responds  to 
inputs  with  a  given  set  of  system  parameters.  These 
parameter  studies  are  used  in  design  and  cost  trade-offs 
which  will  eventually  lead  to  the  optimum  system  design. 
One  way  to  run  parameter  studies  is  to  use  scale  models  and 
testing  chambers.  These  are  primarily  used  in  the  final 
design  stages.  In  preliminary  design  where  there  are  a  lot 
of  variables  computer  simulations  are  best  performed. 

The  computer  model  described  in  Chapter  III  was  used  to 
analyze  some  parameter  changes  that  may  be  associated  with 
ORION.  It  is  expected  that  this  model  will  be  used  to 
determine  the  optimum  configuration  in  terms  of  moments  and 
magnetometer  deployment.  This  chapter  analyzes  the  changes 
in  inertia  ratio  and  spin  rates  as  a  result  of  a  typical 
mission.  It  also  analyzes  the  environmental  effects  of 
gravity  gradient  and  aerodynamic  drag. 

B.  PARAMETER  STUDIES 

The  first  study  was  to  analyze  how  a  change  in  the 
moments  of  inertia  would  affect  the  spin  rate  and  stability. 
The  current  mission  requires  the  deployment  of  four  0.907  kg 
(2  lb)  magnetometers  on  lightweight  booms.   The  booms  were 
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assumed  to  deploy  at  a  constant  rate  of  0.0169  m/sec  to  a 
maximum  of  3.048  m  (10  ft).  The  mass  of  the  booms  was 
considered  to  be  negligible  and  did  not  contribute  to  the 
inertia  of  the  satellite.  The  magnetometers  are  located  in 
the  x-z  plane  of  the  satellite  as  shown  in  Figure  4.1. 

The  effect  of  deploying  the  magnetometers  is  to  increase 
both  Ia  and  If  The  new  moments  of  inertia  can  be 
calculated  by  finding  the  moment  of  inertia  for  the 
satellite  without  the  magnetometers  and  then  adding  the 
moment  of  inertia  for  each  magnetometer.   Thus 

Ia   =   0.5msr2  +  4mm(lb  +  r) 2  (4.1) 

It   =   0.083ms(3r2+h2)  +  2mm(lb+r)2  (4.2) 

where : 

ms   =  mass  of  satellite  minus  the  mass  of  4 
magnetometers 

mm  =  mass  of  one  magnetometer 

r   =   radius  of  the  satellite 

h   =  height  of  the  satellite 

lb   =   distance  of  the  magnetometer  from  the  outside 
of  the  satellite. 

Substituting  in  the  values  for  ORION,  the  moments  of  inertia 
can  be  found  at  any  boom  length  by 

Ia   =   3.203  +  3.628(lb  +  0.242)2  (4.3) 

It   =   8.81  +  1.81(lb  +  0.242)2  (4.4) 
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Figure  4 . 1   Magnetometer  Deployment  Diagram 
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Figure  4.2  shows  this  relationship  between  the  boom  length 
and  the  moments  of  inertia.  Ia  becomes  larger  than  It  at  a 
boom  length  of  1.53  m.  If  energy  could  be  dissipated,  then 
with  Ia  larger  than  It  the  satellite  would  become  stable 
about  the  spin  axis  [Ref.  4:p.  117]. 

Deployment  of  magnetometers  also  has  an  effect  on  the 
spin  rate  and  could  be  used  to  change  it  without  using  the 
spin-up  thrusters.  Using  the  principle  of  the  conservation 
of  angular  momentum,  if  oo^  is  assumed  to  be  zero  then 

H   =   Iau)y  (4.5) 

and  as  Ia  increases  co  „  must  decrease.  Figure  4.3  shows  how 
the  spin  rate  changes  with  boom  length. 

As  noted  earlier,  the  magnetometer  deployment  can  be 
used  to  change  the  spin  rate  of  the  satellite.  If  the  sat- 
ellite is  spinning  at  60  rpm  (wy  =  -6.283  rad/sec)  ,  a  boom 
length  of  1.097  m  will  slow  the  spin  to  20  rpm  (coy  =  -2.094 
rad/sec).  A  1.87  m  deployment  will  slow  the  spin  to  10  rpm 
(coy  =  -1.047  rad/sec).  Adjusting  the  boom  length  to  reach  a 
specific  spin  rate  is  applicable  if  the  mission  does  not  re- 
quire the  magnetometers  a  specific  distance  from  the  satel- 
lite. It  is  more  probable  that  the  boom  must  be  a  specific 
length  and  the  satellite  must  spin  at  a  specific  rate. 

There  are  two  ways  to  achieve  this  configuration.  The 
first  method  is  to  deploy  the  magnetometer  to  the  desired 
length,  then   adjust   the   spin   as   necessary  with   the 
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Figure  4.2   Change  in  Moments  of  Inertia  with  Boom  Length 
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Figure  4.3   Change  in  Spin  Rate  with  Magnetometer  Deployment 
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thrusters.  A  second  method  is  to  despin  the  satellite  with 
the  thrusters  then  deploy  the  magnetometers  and  readjust  the 
spin  with  the  thrusters.  If  the  change  in  spin  and  boom 
length  is  small,  it  may  be  possible  to  despin  the  satellite 
to  a  rate  such  that  the  subsequent  deployment  of  the  booms 
will  despin  the  satellite  to  the  required  rate. 

Since  the  firing  interval  of  the  thrusters  is 
proportional  to  the  amount  of  fuel  used,  it  can  be  used  to 
compare  the  two  methods.  To  illustrate  this,  assume  the 
satellite  is  initially  spinning  at  60  rpm  with  Ia  =  3.239 
kg-m2 .  A  boom  length  of  2  m  and  a  spin  rate  of  20  rpm  is 
required.  At  2  m,  Ia  is  21.447  kg-m2.  Using  the  first 
method,  after  deployment  the  spin  rate  is 


w 


*2  ~  Xa2 

(3.239)  (-6.283) 
21.477 

=   -0.9476  rad/sec  (4.6) 

Assuming  the  total  torque  of  the  thrusters  is  Fxr  =  0.2168 
n-m,  the  firing  interval  dt  is 


dt   = 


I  0  (2.094  -U)  9) 
5l  z     y  *- 

F  xR 


(21.447)  (2.094-0.9476) 
0.2168 

=   113.3  sec  (4.7) 
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Using  the  second  method,  the  firing  interval  reguired  to 
despin  the  satellite  to  20  rpm  is 

dt      (3.239)  (6.283-2.093) 


(0.2168) 
=   62.6  sec  (4-8) 

Deploying  the  magnetometer  to  2  meters  causes  the  spin  rate 
to  decrease  to  0.3152  rad/sec.  The  firing  interval  reguired 
to  readjust  the  spin  is  176.1  sec.  The  total  firing 
interval  for  the  second  method  is  2  3  8.7  sec,  more  than  twice 
that  of  the  first  method. 

This  indicates  that  the  first  method  is  more  fuel 
efficient  in  deploying  the  magnetometers  and  achieving  a 
specific  spin  rate.  A  disadvantage  of  both  these  methods  is 
that  if  the  booms  are  too  flexible,  adjusting  the  spin  may 
cause  the  magnetometers  to  lag  the  boom-satellite  attachment 
points  and  cause  too  great  a  stress  in  the  booms. 

1.   Stability 

Whether  a  is  changed  by  deploying  magnetometers  or 
by  changing  the  internal  mass  distribution  of  the  satellite, 
it  will  have  an  effect  on  the  stability  of  the  satellite. 
Figure  4.4  shows  how  a  change  in  a  affects  A  and  wn.  Since 
the  roots  of  the  open  loop  transfer  function  depend  on  A  , 
they  also  change  with  a  as  was  shown  in  Table  3.1.  Figure 
4.5  shows  the  responses  of  the  satellite  to  a  1.0  n-m  pulse 
of  0.01  sec  as  a  changes.  This  was  done  at  three  different 
spin  rates  of  60  rpm,  2  0  rpm,  and  10  rpm.   The  advantage   of 
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Figure  4.5   Satellite  Response  to  a  Short  Duration  Torque 


59 


increasing  a  is  to  make  the  satellite  less  susceptible  to 
nutation.  Figure  4.6  shows  that  the  higher  the  spin  rate 
and  the  larger  the  inertia  ratio,  the  smaller  the  nutation 
will  be  for  a  given  input  torque. 

C.   ENVIRONMENTAL  EFFECTS 

The  aerodynamic  drag  and  gravity  gradient  torques  were 
modelled  using  the  equations  in  Chapter  III.  These 
equations  were  incorporated  into  the  simulation  to  see  what 
effect  they  would  have  on  the  satellite. 

Since  the  drag  force  is  dependent  on  the  density  of  the 
atmosphere,  it  decreases  as  the  orbit  height  increases. 
This  is  shown  in  Table  4.1.  At  an  altitude  of  200  km  the 
drag  torque  is  1.383  xio"3  n-m  and  decreases  to  1.762  x 10~9 
at  an  altitude  of  1000  km. 

TABLE  4.1 
TYPICAL  VALUES  FOR  AERODYNAMIC  DRAG 

ALT    (km)  Density    (kg/m3)  Drag    (n)  TD    (n-m) 

100                             6    xlO"7  15.84  1.204 

200                             7    xlO"10  .0182  1.383    xio"3 

400                             1   xlO'11  2.523  xlO"4  1.917x  10"5 

1000                              lx    10~15  2.318    xlO"8  1.762   x 10"9 

Gravity  gradient  torque  is  also  dependent  on  orbital 
altitude      since      it      depends   on   co0.      Unlike      drag,      however, 


60 


0.30 


0.40 


r 
1.06 


040  0.7S  0.90 

INERTIA  RATIO  (SIGMA) 


lJW 


1.00 


Figure  4.6   Nutation  Angles  for  Different  Spins  and  Inertias 
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gravity  gradient  torque  is  not  a  constant  but  depends  on  the 
moments  of  inertia  and  the  orientation  of  the  satellite.  At 
an  altitude  of  200  km  the  maximum  torque  is  3.99  x10~9  n-m. 
This  decreases  to  3.72 x  10~15  n-m  at  1000  km. 

An  orbital  altitude  of  2  00  km  was  used  to  determine  the 
effect  environmental  forces  have  on  ORION.  Inertia  ratios 
of  0.364  and  1.488  and  spin  rates  of  10,  20  and  60  rpm  were 
used  to  give  the  widest  ranges  of  results.  These  results 
are  shown  in  Figures  4.7,  4.8  and  4.9.  At  10  rpm  the 
maximum  nutation  is  2.053 x  10~3  radians  for  a  a  of  0.364  and 
5.33  x  10~5  radians  for  a  a  of  1.488.  At  20  rpm  the  nutation 
angles  are  5.173 x  10~4  radians  and  1.17  xio~5  radians  and 
for  60  rpm  they  are  6.0 x  io~5  radians  and  1.33  xlO"6 
radians. 

As  can  be  seen  in  the  figures,  the  nutation  is 
oscillatory.  This  is  because  of  the  oscillatory  nature  of 
the  environmental  torques,  principally  the  gravity  gradient. 
When  the  satellite  nutates  to  its  maximum  angle,  the  torque 
effect  is  reversed  causing  the  nutation  to  decrease.  When 
the  nutation  is  zero,  the  torque  effect  is  again  reversed 
and  the  process  repeats  itself.  The  frequency  of  the  effect 
is  related  to  the  spin  rate  of  the  satellite.  This  is  seen 
in  the  0  term  of  the  torque  equations. 

If  the  accuracy  requirement  of  the  attitude  control 
system  is  small  enough,  environmental  effects  must  be  taken 
into  account  in  the  simulations.   They  will   be   ignored   in 
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Figure  4.8   Nutation  Due  to  Environmental  Torques;  2  0  Rpm 
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further  simulations  described  in  this  thesis  since  their 
effect  is  less  than  the  2  degree  accuracy  requirement  of 
ORION.  This  will  simplify  the  model  without  a  large  error 
in  accuracy. 

D.   ACQUISITION  PHASE  PROFILE 

After  analyzing  the  effects  of  a  changes  and 
environmental  torques,  the  next  step  was  to  simulate  the 
motion  of  the  satellite  from  launch  to  spin-up.  Drag  and 
gravity  gradients  effects  were  included  in  this  simulation. 

The  flight  profile  consisted  of  launching  the  satellite 
with  no  spin.  All  initial  values  of  <J>  ,  ¥  and  0  were  assumed 
to  be  zero.  At  50  seconds  a  step  input  of  Ty  =  0.2168  n-m 
was  applied.  This  caused  the  value  of  u)v  to  increase.  At 
144  seconds  the  input  was  removed.  The  value  of  coy  was 
-6.29  rad/sec  (60  rpm)  .  Figure  4.10  shows  how  a  and  ¥ 
change  during  this  portion  of  the  mission. 

During  the  drift  period  environmental  torques  cause  the 
satellite  to  rotate  0.013  degrees  in  roll  and  11.04  degrees 
in  yaw.  As  the  spin  increases  the  a  and  ¥  begin  an 
oscillatory  motion.  The  amplitude  of  the  motion  depends  on 
the  spin  rate.  As  the  spin  rate  increases,  the  satellite 
becomes  more  resistive  to  torque  effects.  When  the  spin 
rate  becomes  a  constant,  the  amplitude  of  oscillation  due  to 
environmental  torques  becomes  a  constant  although  the  scale 
in  Figure  4.10  is  too  large  to  show  this.  At  144  seconds,  (f> 
is   3.89   degrees   and  ¥  is   12.54   degrees.   The    nutation 
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Figure  4.10   Acquisition  Phase  Roll  and  Yaw 
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caused  by  environmental  torques  is  0.0015  degrees,  much  less 
than  the  acceptable  nutation  of  2  degrees.  After  400 
seconds  with  no  controlling  inputs,  $  is  7.4  degrees,  ¥  is 
10.8  degrees  and  the  nutation  is  0.0015  degrees. 

From  this  point,  the  magnetometers  could  be  deployed  or 
control  torques  could  be  applied  to  reorient  the  satellite. 
Chapter  V  discusses  the  methods  of  controlling  nutation  and 
reorienting  the  satellite. 
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V.   MISSION  PHASE  SIMULATION 

A.   SATELLITE  CONTROL 

The  process  of  maintaining  the  correct  orientation  of  a 
satellite  consists  of  eliminating  the  precession  and 
nutation  effects  caused  by  external  or  internal 
disturbances.  The  simplest  case  is  when  the  spin  axis  and 
angular  momentum  vector  are  aligned  and  are  perpendicular  to 
the  orbital  plane. 

This  presents  a  two  part  control  problem.    The  first 

part  is  to  reduce  nutation  by  forcing  cot  equal  to  zero.   The 

second  part  is  maintaining  the  angle  between  the  satellite's 

angular  momentum  vector  H  and  the  orbital  angular  momentum 

vector  H0.   The  relationship  between  H  and  H0  is 

H  -H 
o 
Cos  a  - 


|H0||H| 

The  angle  a  can  be  obtained  by  precessing  the  satellite  with 
control  torques  and  then  eliminating  any  resulting  oot. 

co  t  can  be  forced  to  zero  by  applying  a  control  torque  of 
magnitude  equal  to  Itwt  ^n  the  direction  of  Itwf  wt  can  be 
calculated  from  measurements  of  $,  V  ,  4>  and  V  .  The 
magnitude  of  torque  required  is  determined  by 

F  *  r  *  dt   =   Itwt 
where: 
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F  =   the  impulsive  force  of  the  thruster 

r  =  distance  of  the  thruster  from  the  center 
of  mass 

dt   =   time  over  which  the  thruster  is  fired. 


An  ideal  thruster  would  provide  an  impulsive  force  at  the 
correct  time  in  the  spin  cycle  and  dt  would  be  zero.  The 
firing  profile  of  an  actual  thruster  can  be  approximated  by 
the  trapezoid  model  shown  in  Figure  5.1  [Ref.  3: p.  21].  t-^ 
is  the  time  the  thruster  is  commanded  on,  t2  is  the  time  the 
thruster  reaches  maximum  force,  t4  is  the  time  the  thruster 
is  commanded  off  and  t5  is  the  time  the  thruster  force  is 
zero.  The  total  time  dt  is  from  t^  to  t$.  tc  is  the  time 
at  which  the  effective  torque  is  applied. 

t2  and  t5  can  be  determined  from  experimental  or  the 
manufacturer's  data.  Knowing  tc,  t2  and  t5,  t±  and  t4  can 
be  determined  by 

Tan  Atc  =  b/a  (5.1) 

(jO  I  /  c  2  \ 

-^  =  a  cos  At  +  b  sin  At  {         } 

r  c         c 


where: 


cos  At2  -  cos  At,  cos  At  -  cos  At4 
a   =      A(t2-tx)  A(t5-t4) 

sin  At2  -  sin  At,  sin  At^  -  sin  At4 
b   "      ACt^)  A(t5-t4) 
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Figure  5.1   Typical  Thruster  Profile  (Trapezoid  Model) 
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Figure  5.2  shows  the  relationships  between  the  firing 
profile  and  timing. 

An  even  simpler  approximation  is  the  pulse  function. 
This  would  occur  if  t2  =  t±  and  t4  =  t5.  Knowing  that  dt  = 
t4-t2,  tc  =  (l/2)dt.   Therefore 

t2   =   tc  -  (1/2)  dt  (5.3) 

t4   =   tc  +  (l/2)dt  (5.4) 

Once  the  magnitude  of  control  torgue  and  the  firing  arc 
are  known,  the  problem  is  when  during  the  spin  cycle  to  fire 
the  thruster  to  eliminate  ^^.  The  acceleration  of  any  point 
(x,y,z)  on  the  satellite  is 

2 
ax   =   -xo)s  +  wnyA  sin  At  (5.5) 

ay   =   (<%"X)A(z  cos  At  +  x  sin  At)  (5.6) 

2 
az   =   -Zcos  +  conyAcos  At  (5.7) 

where  A  is  proportional  to  the  nutation  angle  [Ref.  4: p. 
116].  If  an  accelerometer  is  placed  in  the  x-z  plane,  the 
axial  acceleration  reaches  a  positive  peak  when  ^t  points 
towards  the  accelerometer.  Placing  a  thruster  90  degrees 
from  the  accelerometer  produces  a  torgue  opposite  to  co^.  If 
the  accelerometer  is  located  on  the  z  axis, 

ay   =   (ws"A)Az  cos  At  (5.8) 

and   the   peak   acceleration   is   reached  when   day/dt  =  0. 

72 


Figure  5.2   Relationship  Between  Spin  Rate  and 
Firing  Times 
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Since 


dav/dt   =   -(ws-A)AzA  sin  At  (5.9) 


then 


day/dt      =      0      when      t      =      tt/A 

The  value  of  acceleration  at  the  firing  point  is 

aymax   =   -(ws-A)zA  (5.10) 

Using  the  value  of  t  =  tt/A  as  tc,  the  on  and  off  times  for 
the  thruster  can  be  calculated.  If  the  spin  rate  is  too 
fast,  dt  may  encompass  several  revolutions  of  the  satellite. 
To  provide  a  torque  in  a  given  direction,  however,  the 
firing  must  be  limited  to  an  arc  centered  about  tc. 
Therefore,  the  total  firing  time  required  (dt)  must  be 
divided  into  segments  occurring  at  integer  multiples  of  tc. 

For  example,  assume  the  nutation  angle  is  1.0  deg 
(0.0174  rad)  ,  It  is  8.91  kg-m2 ,  coy  i-s  -6.283  rad/sec,  and  H 
is  2  0.35  kg-m2/sec2.   Therefore 

Itoot  =  H  sin  n 

=   0.355  kg-m2/sec2  (5.11) 

Assuming  F  is  0.448  n  and  R  is  0.242  m 

It03t 
dt       =       FT 

=      3.28    sec    .  (5.12) 
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At  60  rpm  the  satellite  makes  one  revolution  per  second.  If 
the  firing  arc  is  90  degrees,  the  firing  time  per  arc  is 

2^-   =   0.25  sec  (5.13) 

Y 

Therefore,  13  firings  will  be  reguired  to  eliminate  the 
nutation.  If  the  spin  rate  was  10  rpm,  only  2  firings  would 
be  reguired.  Each  firing  must  occur  at  aymax  to  avoid 
adding  energy  to  the  system. 

This  method  of  active  nutation  control  was  simulated 
using  the  computer  model.  Nutation  was  started  by  inputting 
a  5  n-m  pulse  of  0.01  seconds.  After  0.03  seconds  the 
control  torgue  could  be  applied.  The  minimum  reguired 
nutation  was  set  at  0.001  radians  to  avoid  computational 
errors.   The  firing  arc  was  limited  to  45  degrees. 

To  apply  the  correct  control  torgue,  the  program 
calculated  the  value  of  ay  which  would  cause  the 
acceleration  peak  to  be  in  the  center  of  the  firing 
interval.  When  the  calculated  value  of  ay  was  greater  than 
this  value,  Tz  was  set  egual  to  2Fr/dt  where  dt  is  the 
firing  arc  time  period.  As  long  as  the  nutation  was  greater 
than  0.001,  T2  was  applied  at  the  peak  positive  acceleration 
point. 

This  simulation  was  run  for  3  different  spin  rates  of 
10,  20  and  60  rpm  at  an  inertia  ratio  of  0.363.  It  was  also 
run  for  two  different  inertia  ratios  and  a  constant  spin 
rate.    These  results  are  shown  in  3  series  of  3  figures. 
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Each  series  shows  the  nutation  angle  and  applied  control 
torque,  the  $  and  ^  angles  and  the  firing  diagrams  for  a 
different  spin  rate. 

In  Figure  5.3,  the  spin  rate  was  10  rpm.  The  initial 
nutation  was  1.45  *  10~3  radians.  It  was  reduced  to  0.001 
radians  at  4.4  seconds.  Tz  was  applied  to  reduce  the 
nutation  to  0.001  radians.  From  Figure  5.4,  cj)  and  y  were 
reduced  to  steady-state  values  of  1.67  x 10~2  radians  and 
-1.66  x  10~4  radians.  The  values  of  <$>  and  ¥  oscillated  about 
these  values  because  the  nutation  was  not  completely 
eliminated.  Figure  5.5  shows  the  firing  diagram  for  this 
case.   Tz  was  8.3  3 x  10~2  n-m  over  a  firing  time  of  .6  sec. 

In  Figure  5.6,  the  spin  rate  was  changed  to  2  0  rpm  which 
resulted  in  an  initial  nutation  of  7.28 x  10~3.  The  first  of 
3  control  firings  occurred  at  2.2  seconds.  Nutation  was 
reduced  after  12  seconds.  In  Figure  5.7,  the  values  of  (J) 
and  y  are  7.91  x  10~3  radians  and  -2.1  ><10~4  radians.  From 
the  firing  diagram  of  Figure  5.8,  Tz  was  4.1 x  10~2  n-m  over 
a  pulse  of  .36  sec.  Three  pulses  were  needed  for  a  total 
firing  time  of  1.08  seconds. 

In  Figure  5.9,  a  spin  rate  of  60  rpm  was  used.  The 
initial  nutation  was  2.43 x  10~3  radians,  the  first 
application  of  control  occurred  at  .7  seconds  and  the 
nutation  was  reduced  after  27.5  seconds.  <jj  and  T  from 
Figure  5.10  were  2.27  x  10~3  and  -2.0  x  10~5  radians.  From 
the  firing  diagram  of  Figure  5.11,  18  applications   of  Tz 
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Figure  5.4   Roll  and  Yaw  Angles;  wy  =  10  Rpm 
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Figure  5.5   Firing  Diagram;  wy  =  io  Rpm 
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Figure  5.6   Nutation  Control;  coy  =  20  Rpm 
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Figure  5.8   Firing  Diagram;  iOy  =  20  Rpm 
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were  required.  Tz  was  1.3  xio"2  n-m  with  a  pulse  width  of 
.23  seconds.  A  total  firing  time  of  4.2  seconds  was 
required. 

Figure  5.12  shows  the  difference  in  control  for  the  same 
spin  rate  of  2  0  rpm  but  for  two  different  inertia  ratios. 
When  o  =  0.363  it  required  3  applications  to  reduce  the 
nutation.  When  was  increased  to  1.12  the  initial  amount 
of  nutation  was  less  and  only  one  application  of  control 
torque  was  required. 

As  a  result  of  these  simulations,  it  can  be  seen  that 
the  slower  spin  rates  require  less  firing  times  and  number 
of  pulses.  This  is  because  the  firing  arc  time  is  larger. 
A  low  number  of  pulses  also  reduces  the  probability  of  a 
timing  error  which  may  cause  nutation  instead  of  eliminate 
it.  In  every  case  shown,  <j>  reached  a  steady  state  value 
while  y  approached  zero.  This  is  a  result  of  using  Tz  for 
the  input  and  control  torques. 

Fewer  pulses  resulted  in  shorter  firing  times  and 
quicker  responses.  This  equates  to  a  fuel  savings.  The 
model  can  also  be  run  with  different  size  thrusters  by 
changing  the  value  of  F  in  the  program.  Larger  firing  arcs 
would  also  decrease  the  response  times  especially  at  higher 
spin  rates. 

B.   SATELLITE  REORIENTATION 

The  control  system  may  be  used  to  reorient  the  satellite 
as   well    as   reduce  the   nutation.   Reorientation   can 
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Figure  5.12   Nutation  Variation  with 3  ; toy  =  20  Rpm 
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be  accomplished  either  as  part  of  the  nutation  reduction  or 
as  a  separate  control  sequence.  It  was  shown  in  Figures 
5.4,  5.7  and  5.10  that  although  the  nutation  was  reduced,  <!> 
and  4/  were  not  returned  to  their  initial  values. 

Reorientation  during  nutation  control  is  a  two-step 
process.  The  first  step  involves  reducing  the  nutation 
angle  by  half  its  original  amount.  This  will  cause  a  new 
nutation  cone  to  be  established.  When  the  spin  axis  moves 
180°  along  the  new  cone,  control  torques  are  applied  to 
reduce  the  nutation  to  the  minimum  level. 

Because  of  the  differences  between  nutation  frequency, 
spin  rate  and  A  ,  the  thrusters  may  not  be  aligned  with  the 
cot  vector  when  the  satellite  reaches  the  180°  point  on  the 
nutation  cone.  The  time  required  for  the  spin  axis  to 
travel  180°  is 

At   =   7T/con  (5.14) 

The  satellite  will  rotate  about  its  axis  through  an  angle 

a  =   wyAt  (5.15) 

in  this  time  period.  The  relative  rotation  of  the  oot  vector 
is 

3   =   A At  (5.16) 

over  the  same  time  period.  The  two  angles  a  and  3  are  equal 
only  if  a  =  1.0. 
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Figure  5.13  shows  the  relationship  between  cuy,  wR  and  X. 
Since 

X      =      (a-l)coy 

=   (Ia/It  "  !)^y  (5.17) 

and 

wn      =      Ia/Itajy  (5.18) 

The   two   angles    a  and    3  can  be   calculated   as 

3       =      t*  (Q~1^  (5.19) 

a      =      tt/o  (5.20) 

Since  over  the  required  time  period  the  thrust  vector  may 
not  be  aligned  with  the  wt  vector,  the  application  of  a 
control  torque  would  add  to  the  nutation. 

To  overcome  this  difficulty,  four  sets  of  thrusters  can 
be  used.  If  the  thrusters  are  located  on  the  x  and  z  axes, 
a  combination  of  thrusters  can  be  fired  to  make  the 
resulting  torque  collinear  with  the  cjo-j.  vector. 

This  is  accomplished  in  the  program  by  saving  the  time 
at  which  the  nutation  is  halved  as  FTIME.  When  the  program 
time  reaches  FTIME  plus  At,  6  is  calculated  and  the  correct 
combination  of  Tx  and  Tz  used  to  reduce  the  nutation  to  an 
acceptable  level  and  restore  cf>  and  ¥  to  their  initial 
values. 
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The  second  method  is  to  reorient  the  satellite  assuming 
there  is  no  nutation.  The  degree  of  reorientation  was  given 
earlier  in  the  chapter  as 

H  -H 

cos  a  =  jtrmr 

Knowing  the  angle  a  and  H0,  the  H  reguired  to  make  a  zero 
can  be  calculated.  This  gives  the  reguired  co-j-.  A  control 
torgue  egual  to  Itco  t  is  applied  to  start  the  satellite 
nutating.  When  A  t  is  egual  to  TT/con  a  control  torgue  egual 
and  opposite  to  I-t^t  ^s  applied  to  reduce  the  nutation. 
This  is  done  in  a  manner  similar  to  that  of  reorienting  the 
satellite  with  nutation. 
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VI.   CONCLUSIONS 

This  thesis  derived  the  equations  which  describe  the 
rotational  motion  of  a  satellite.  It  also  derived  the 
equations  describinq  the  effects  of  qravity  qradient  torque 
and  aerodynamic  draq.  These  equations  were  then  used  to 
construct  a  computer  simulation  model.  The  model  was  used 
to  analyze  how  chanqes  in  moments  of  inertia  and  spin  rates 
affect  the  motion  of  the  satellite. 

It  was  shown  that  extendinq  maqnetometers  on  booms  would 
cause  Ia  and  I-j-  to  increase.  This  chanqed  the  inertia  ratio 
a.  The  spin  rate  of  the  satellite  could  then  be  chanqed  by 
chanqinq  the  boom  lenqth  of  the  maqnetometer.  As  a 
increased,  the  stability  of  the  satellite  was  chanqed. 
A  hiqh  value  of  a  made  the  satellite  resistive  to  nutation 
by  external  forces. 

The  effects  of  qravity  qradient  and  aerodynamic  draq 
were  analyzed.  The  amount  of  nutation  caused  by  environment 
depended  on  the  spin  rate  and  a .  The  maximum  amount  of 
nutation  occurs  at  a  low  spin  rate  and  low  a .  The  nutation 
is  on  the  order  of  2.053  x 10"3  radians  for  a  a  of  0.363  and 
a  10  rpm  spin  rate.  The  effects  of  the  environment  need 
only  to  be  taken  into  account  in  the  motion  analysis  if  the 
accuracy  required  by  the  control  system  is  small.  The 
accuracy  for  ORION  is  to  be  within  2  deqrees  so  qravity 
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gradient  and  drag  were  only  applied  during  the  example 
flight  profile. 

The  control  system  that  is  used  is  an  active  nutation 
control.  An  example  of  how  the  control  system  would  be  used 
was  given.  Assuming  a  firing  arc  of  45°,  4.4  seconds  were 
required  to  reduce  the  nutation  at  a  spin  rate  of  10  rpm,  12 
seconds  for  a  spin  rate  of  20  rpm  and  27.5  seconds  for  a 
spin  rate  of  60  rpm.  The  procedure  for  reorienting  the 
satellite  to  any  desired  attitude  was  also  described. 

Although  the  computer  model  was  written  to  simulate  the 
satellite  as  closely  as  possible,  many  areas  need  more 
detailed  modelling.  The  satellite  was  assumed  to  be  a 
cylinder  of  uniform  mass.  The  actual  mass  distribution 
needs  to  be  determined  and  the  moments  of  inertia 
calculated.  This  would  change  the  inertia  ratio  and 
possibly  cause  some  cross-products  of  inertia  to  develop. 

Another  factor  that  was  not  included  in  this  model  is 
the  effect  of  liquids  slosh  in  the  fuel  tanks.  This  could 
cause  internal  torques  that  adversely  affect  the  motion. 
These  tanks  could  also  be  designed  to  dissipate  energy.  If 
this  could  be  done,  then  when  a  was  greater  than  one,  the 
satellite  could  be  made  stable. 

The  firing  profile  of  the  control  thruster  was  assumed 
in  the  model  to  be  a  pulse  function.  Entering  the  profile 
as  a  trapezoid  would  increase  the  accuracy  of  the  model. 
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When  real  thruster  data  becomes  available  it  could  then  be 
entered. 

As  has  been  stated  before,  many  of  the  variables  in  the 
program  can  easily  be  changed  in  parameter  runs.  This 
allows  a  great  deal  of  flexibility  in  analyzing  design 
changes  as  they  occur. 
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APPENDIX 
ATTITUDE  ANGLE  CONTROL  PROGRAM 


TITLE  ATTITUDE 

*  THIS  PROGRAM 

*  ORION  SATELLI 

*  SATELLITE  MOT 

*  FREQUENCIES  A 

*  THE  ENVIRONME 

*  COMPUTES  THE 

*  NECESSARY  FOR 

*  PROGRAM  CAN  S 


ANGLE  CONTROL  PROGRAM 

PROVIDES  THE  SIMULATION  MODEL  FOR  THE  CONTROL  OF  THE 
TE.  IT  USES  THE  EOUATIONS  DERIVED  IN  CHAPTER  3  TO  MODEL 
ION  AND  ENVIRONMENTAL  TORQUES.   IT  CALCULATES  THE  ANGULAR 
ND  THE  ATTITUDE  ANGLES.   IT  CALCULATES  THE  MAGNITUDE  OF 
NTAL  TOROUES  BASED  ON  0RI3TAL  PARAMETERS.  IT  ALSO 
NUTATION;  ANGULAR  MOMENTUM  VECTOR  AND  THE  CONTROL  TORQUES 

THE  ACTIVE  NUTATION  CONTROL  SYSTEM.   CHANGES  TO  THE 
n  MADE  VIA  PARAM  STATEMENTS  OR  DIRECTLY  IN  THE  LISTING. 


*  VARIABLE  DESCRIPTION 


=  DISTANCE  FROM  THE  CENTER  OF  PRESSURE  TO  THE  CENTER  OF  MASS 

=  THRUST  FORCE  OF  ONE  CONTROL  THRUSTER 

=  RADIUS  OF  THE  SATELLITE 

=  TRANSVERSE  MOMENT  OF  INERTIA 

■  MOMENT  OF  INERTIA  ALONG  THE  SPIN  AXIS 

■  ORBITAL  ALTITUDE  OF  THE  SATELLITE 
=  RADIUS  OF  THE  EARTH 

=  DENSITY  OF  AIR  AT  THE  ORBITAL  ALTITUDE 

■  ANGULAR  VELOCITY  OF  THE  SATELLITE  WRT  EARTH 

■  DRAG  FORCE 

'  ANGULAR  VELOCITY  ABOUT  THE  X  AXIS  WRT  EARTH 

■  ANGULAR  VELOCITY  ABOUT  THE  Z  AXIS  WRT  EARTH 

■  ANGULAR  VELOCITY  ABOUT  THE  Y  AXIS  WRT  EARTH 
=  SATELLITE  ROTATION  ABOUT  THE  X  AXIS 

=  SATELLITE  ROTATION  ABOUT  THE  Z  AXIS 

:  SATELLITE  ROTATION  ABOUT  THE  Y  AXIS 

:  ROLL  ANGLE 

■  YAW  ANGLE 

:  PITCH  ANGLE 

■  DRAG  TORQUE  IN  THE  X  DIRECTION 

■  DRAG  TORQUE  IN  THE  Z  DIRECTION 
=  TOTAL  DRAG  TOROUE 

=  GRAVITY  GRADIENT  IN  THE  X  DIRECTION 

■  GRAVITY  GRADIENT  IN  THE  Z  DIRECTION 
:  TOTAL  GRAVITY  GRADIENT  TORQUE 

■  CONTROL  TOROUE  IN  THE  X  DIRECTION 
'  CONTROL  TORQUE  IN  THE  Z  DIRECTION 
=  CONTROL  TORQUE  IN  THE  Y  DIRECTION 
=  TRANSVERSE  ANGULAR  VELOCITY 

■  TRANSVERSE  COMPONENT  OF  THE  ANGULAR  MOMENTUM  VECTOR 
'  AXIAL  COMPONENT  OF  THE  ANGULAR  MOMENTUM  VECTOR 

■  MAGNITUDE  OF  THE  ANGULAR  MOMENTUM  VECTOR 
=  NUTATION  ANGLE 

'  INERTIA  RATIO 

=  RELATIVE  VELOCITY  OF  WT  WRT  THE  SATELLITE  ROTATION 

'  NUTATION  FREQUENCY  WRT  INERTIAL  SPACE 

=  AXIAL  ACCELERATION  OF  THE  SATELLITE  DURING  NUTATION 

=  FIRING  ARC  OF  THE  CONTROL  THRUSTERS 

•  TIME  OF  PEAK  POSITIVE  AXIAL  ACCELERATION 

=  ACCELERATION  AT  THE  THRUSTER  TURN  ON  TIME 

■  TIME  OF  THE  THRUSTER  FIRING 

=  TIME  FOR  THE  SPIN  AXIS  TO  TRAVEL  PI  RADIANS 
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ALT  =  200.00 

ROE  =  7.0D-10 

RO  =  (6378.0+ALT)*1.0D03 

WO  =  SQRT(3.987D14/RO**3) 

D  =  1.7D14*R0E/R0 

FLAG  =  0 
INCON  WXO=0 . 0 , WZO=0 . 0 , WY0=-2 . 093 , PHIO=0 . 0 , PSIO=0 . 0 , THO=0 . 0 
DERIV 

WXD  =  (TX+TDX+GX-(IT-IA)*WY*WZ)/IT 

WZD  =  (TZ+TDZ+GZ-(IA-IT)*WY*WX)/IT 

WYD  =  TY/IA 

PHID  =  WX  +  PSI*WO 

PSID  =  WZ-PHI*WO 

THD  =  WY+WO 

WX  =  IMTGRL(WXO,WXD) 

WZ  =  INTGRL(WZO,WZD) 

WY  =  INTGRL(WYO,WYD) 

PHI  =  INTGRL(PHIO,?HID) 

PSI  =  INTGRL(PSIO,PSID) 

TH  =  INTGRL(THO,THD) 
DYNAMIC 

*  ENVIRONMENTAL  TOROUE  CALCULATION;  USE  WHEN  DISTRUBANCE  TORQUES  ARE 

*  REQUIRED 

*  TDX  =  -D*L*SIN(TH)*C0S(PSI) 

*  TDZ  =  D*L*COS(TH)*COS(PSI) 

*  TD  =  (TDX**2+TDZ**2)**0.5 

*  GX  =  3.0*(WO**2)*(IT-IA)*SIN(PHI)*COS(PHI)*(COS(TH))**2 

*  GZ  =  3.0*(WO**2)*(IT-IA)*SIN(TH)*COS(TH)*SIN(PHI) 

*  GG  =  (GX**2+GZ**2)**0.5 

k  k  k  k  k  k  k 

*  USE  THE  NEXT  4  STATEMENTS  WHEN  NO  DISTURBANCE  TORQUES  ARE  REQUIRED 

TDX  =  0.0 
TDZ  =0.0 
GX  =  0.0 
GZ  =0.0 

kkkkkkk 

TX  =  0.0 

TY  =  0.0 
* 

*  STEP  FUNCTION  USED  TO  SET  AN  INITIAL  NUTATION  VALUE 

TZ  =  5.0*STEP(0.01)-5.0*STEP(0.02) 

k 

*  CALCULATION  OF  ANGULAR  VELOCITIES,  MOMENTUM  VECTORS  AND  FREQUENCIES 

WT  =  (WX**2  +  WZ**2)**0.5 

HT  =  IT*WT 

HS  =  IA*WY 

H  =  (HT**2  +  HS**2)**0.5 

NU  =  ARSIN(HT/H) 

SIGMA  =  IA/IT 

LAMDA  =  (SIGMA-1.0)*WY 

WN  =  SIGMA*WY 

k 

*  CALCULATION  OF  THE  FIRING  TIME  AND  FIRING  ARC 

AY=  (WY-LAMDA)  *  R  *  COS (LAMDA*TIME) 

FARC  =  0.392699/ABS(WY) 

TC  =  3.14159/LAMDA 

AON  =  (WY-LAMDA)*R*COS(LAHDA*(TC-FARC)) 

k 

*  THIS  STATEMENT  ALLOWS  THE  NUTATION  TO  BE  PRESENT  BEFORE  CONTROL  OCCURS 

IF  (TIME. LE. 0.2)  GO  TO  15 

k 

*  CONTROL  CALCULATION  AND  FIRING.  IF  FLAG  IS  0  THE  CONTROL  SYSTEM  WILL 

*  REDUCE  THE  INITIAL  NUTATION  BY  HALF  THEN  FIRE  TO  REDUCE  THE  ATTITUDE 

*  ANGLES  TO  ZERO.  IF  THE  FLAG  IS  1  THE  CONTROL  SYSTEM  WILL  REDUCE  THE 

*  NUTATION  TO  NUREF  ONLY. 

IF ((NU-NUREF).LE. 0.00001)  THEN 

TX  =0.0 

TZ   =   0.00000 
ELSE 

IF(FLAG.EQ.O)THEN 
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END 
IF 


DNU  =  (NU-NUREF)/2.0  +  NUREF 
FLAG  =  1 
END  IF 
I F ( AY . GT . AON . AND . FLAG . EQ . 1 ) THEN 

FT I ME  =  TIME 

TZ  =  2.0*F*R*FARC 

TX  =  0.0 
ELSE 

TZ=0.0 

TX  =0.0 

(NU.LE. DNU) THEM 

CANCEL  =  FTIME  +  ABS (3 . 14159/WN) 
FLAG  =  2 
END  IF 

*  CALCULATION  OF  THE  CORRECT  THRUSTER  COMBINATION  REQUIRED  TO  REDUCE 

*  THE  ATTITUDE  ANGLES  TO  A  MINIMUM. 

IF  (TIME . GE . CANCEL . AND . FLAG . EQ . 2 ) THEN 
DT  =  CANCEL  +  2*FARC 
CONT  =  4*F*R*FARC 

TCL  =  -CONT*STEP( CANCEL)  +  CONT*STEP (DT) 
IF  ( SIGMA. GT. 1.0)  THEN 

BETA  =  ((SIGMA-1.0)*3.14159/SIGMA)  -6.283 
ELSE 

BETA  =  (SIGMA-1.0)*3.14159/SIGMA 
END  IF 
IF  (BETA. GT. 0.0. OR. BETA. LE. -1.5708)  THEN 

DELTA  =  ABS(EETA) 

TZ  =  -TCL*COS (DELTA) 

TX  =  -TCL*SIN(DELTA) 
ELSEIF  (BETA. GT. -1.5708. OR. BETA. LE. -3. 14159)  THEN 

DELTA  =  3.14159-  ABS(BETA) 

TZ  =  TCL*COS (DELTA) 

TX  =  -TCL*SIN(DELTA) 
ELSEIF  (BETA. GT. -3. 14159. OR. BETA. LE. -4. 7123)  THEN 

DELTA  =  ABS(BETA)  -  3.14159 

TZ  =  TCL*C03(DELTA) 

TX  =  TCL*SIN(DELTA) 
ELSEIF  (BETA. GT. -4. 7123. OR. BETA. LE. -6. 283)  THEN 

DELTA  =  6.283  -  ABS(EETA) 

TZ  =  -TCL*COS (DELTA) 

TX  =  TCL*SIN(DELTA) 
END  IF 
IF 


END 
15   END  IF 
COMTRL  FINTIM=20.0,  DELT=0.01 
SAVE  .01,  NU,  PHI,  PSI,  WX,  WZ ,  WY,  SIGMA, 
H,  HS,  HT,  PHID,  PSID,  THD ,  TH, 


LAMDA,WN,AY,AON, 
TX,  TY,  TZ,  WT 


GRAPH  (G1,DE=TEK618)  TIME(UN=SEC)  PHI (UN=RAD ,LO=-3 . OE-03 , SC=1 . 5E-03) 


GRAPH 

GRAPH 

LABEL 

LABEL 

LABEL 

LABEL 

END 

STOP 


PSI(UN=RAD,LO=-3.0E-03,SC=1.5E-03) 
TIME  NU(UN=RADS) 
TIME  TX  TZ(LO=-32E-03,SC=.0165) 
ROLL  AND  YAW  ANGLES;  WY  =  20RPM 
YAW  ANGLE;  SIGMA=0.3635  &  1.479;  WY=-6.283 

G3;G9,DE=TEK613)  NUTATION  ANGLE;  WY  =  20  RPM 

G4,DE=TEK618)  FIRING  DIAGRAM  WY  =  20  RPM 


'G3,DE=TEK618 
G4,DE=TEK618 
G1,DE=TEK618 
G2,DE=TEK618 
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